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ABSTRACT 

Protoplanetary discs are now routinely observed and exoplanets, after the numerous 
indirect discoveries, are starting to be directly imaged. To better understand the planet 
formation process, the next step is the detection of forming planets or of signposts of 
young planets still in their disc, such as gaps. A spectacular example is the ALMA 
science verification image of HL Tau showing numerous gaps and rings in its disc. 

To study the observability of planet gaps, we ran 3D hydrodynamical simulations 
of a gas and dust disc containing a 5 Mj gap-opening planet and characterised the spa¬ 
tial distribution of migrating, growing and fragmenting dust grains. We then computed 
the corresponding synthetic images for ALMA. For a value of the dust fragmentation 
threshold of 15 ms“' for the collisional velocity, we identify for the first time a self- 
induced dust pile up in simulations taking fragmentation into account. This feature, 
in addition to the easily detected planet gap, causes a second apparent gap that could 
be mistaken for the signature of a second planet. It is therefore essential to be cautious 
in the interpretation of gap detections. 
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1 INTRODUCTION 

The holy grail in the field of planetary science is to catch 
planet formation in action, and the most anticipated sign¬ 
post of massive planets in discs are gaps. Many gas-rich pro¬ 
toplanetary discs have now been observed with gaps and 
inner cavities, especially at sub-millimetre and millimetre 
wavelengths (e.g. Casassus et al. 2013; Canovas et al. 2015; 
van der Marel et al. 2015). These so-called transitional discs 
have inner regions which have undergone substantial clear¬ 
ing, due to a combination of grain growth, photoevaporation 
and/or interaction with an embedded planet or companion 
(Espaillat et al. 2014). While there are now a few examples 
of discs which show tantalising evidence of embedded pro¬ 
toplanets and disc clearing (e.g. LkCa 15, Kraus & Ireland 
2012; HD 100546, Quanz et al. 2015) direct observational ev¬ 
idence of planet clearing remains elusive. Planet gaps have 
been predicted theoretically for several decades (Papalolzou 
& Lin 1984) and many numerical studies of planets opening 
gaps have been conducted (e.g. Kley 1999; Paardekooper & 
Mellema 2004; de Val-Borro et al. 2006; Crida et al. 2006; 
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Fouchet et al. 2007; Ayliffe et al. 2012), as have predictions 
of their observability (Wolf et al. 2002; Gonzalez et al. 2012; 
Ruge et al. 2013). 

More recently, the spectacular ALMA science verifica¬ 
tion images of HL Tau (ALMA Partnership et al. 2015) sur¬ 
prised everyone with its clearly discerned multiple rings and 
gaps, which even more surprisingly are axisymmetric. It is 
tempting to attribute each ring plus gap in the ALMA im¬ 
age of HL Tau to a planet, but is this the right thing to do? 
Zhang et al. (2015) invoke potential accumulation of dust at 
condensation fronts. However, the location of the first three 
gaps observed by ALMA in HL Tau are found to be close to 
mean motion resonances, suggesting that the rings could be 
generated by planets (ALMA Partnership et al. 2015). The 
question we address in this Letter is whether all gaps seen in 
protoplanetary discs are always planet gaps, or could they 
result from alternative processes occurring during the early 
stages of planet formation? We explore here the effects of 
dust pile up as such an alternative. 

Millimetre continuum images trace the spatial distribu¬ 
tion of dust grains. Their evolution in discs is governed by 
two parameters. The first is the Stokes number, St, which 
is the ratio of the drag stopping time to the orbital period. 
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Figure 1. Maps of the volume density of the gas (left) and dust (centre) and of the average dust grain size (computed with the SPH 
kernel, right) in the disc midplane after 100 000 yr. 
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The value of St determines the efficiency of vertical settling, 
radial drift and growth/fragmentation of grains in the disc 
(e.g. Weidenschilling 1977; Nakagawa et al. 1986; Stepinski 
& Valageas 1997; Youdin & Shu 2002; Garaud & Lin 2004; 
Haghighipour 2005; Brauer et al. 2008; Birnstiel et al. 2010; 
Laibe 2014). In particular, grain migration is fastest when 
St ~ 1. St is a function of the grain size and density, but 
depends also on the local gas conditions. As such, it varies 
as a function of the distance to the central star for a given 
grain size. HL Tau is a young and massive extended disc 
supplied by in-falling material from the surrounding enve¬ 
lope. For such discs, millimetre-sized grains have St ~ 1 at a 
few tens of au, at which point the dust and gas strongly de¬ 
couple (Laibe et al. 2012). The second parameter involved is 
the dust-to-gas ratio e. For e of order unity, dust affects the 
gas as much as gas affects dust due to aerodynamic drag. 
Although the dust backreaction due to dust inertia is neg¬ 
ligible in most astrophysical systems, it becomes critical in 
regions of discs where dust concentrates, such as particles 
traps. 

In this Letter we present a mechanism where grain 
growth and fragmentation result in dust distributions that 
produce apparent gaps when the dust backreaction is taken 
into consideration. We provide the details on our numerical 
simulations and synthetic observables in Sect. 2, discuss the 
origin of the apparent gap in Sect. 3, and conclude in Sect. 4. 

2 METHODS AND RESULTS 
2.1 Hydrodynamic simulations 

We study the evolution of gas and dust in protoplanet ary 
discs using our 3D, two-phase, Smoothed Particle Hydro¬ 
dynamics (SPH) code, described in Barriere-Fouchet et al. 
(2005). We self-consistently treat the aerodynamical drag 
between gas and dust, including the backreaction of dust on 
gas. We follow the standard SPH implementation of artifi¬ 
cial viscosity (Monaghan 1989) and use cTsph 0.1 and/lspH 0.5, 
corresponding to a uniform Shakura & Sunyaev (1973) pa¬ 
rameter ass ^ 10 ^ (see Fouchet et al. 2007, for a discussion). 
Arena & Gonzalez (2013) showed that the expected prop¬ 
erties of turbulence are naturally reproduced by the SPH 
formalism. Grain growth is treated as described in Laibe 
et al. (2008), via the prescription developed by Stepinski & 
Valageas (1997) based on the determination of the relative 


velocity of grains, V^gi x Cg sJassSt/l St, as a function of lo¬ 
cal properties of the gas phase (e.g. its sound speed Cs) and 
grain size. We implement grain fragmentation by defining 
a velocity threshold, Vfrag, to which we compare Fpel- Golli- 
sions occurring at speeds below the threshold velocity lead 
to grain growth, while collisions at greater speeds result in 
the shattering of grains, which we model by decreasing the 
size of the representative grains. The numerical procedure is 
detailed in Gonzalez et al. (2015). 

We model a disc typical of Glassical T Tauri Stars 
(CTTS) of mass M^^sc 0.01 Mq around a 1 Mq star, with 
an initially uniform dust-to-gas ratio e 10“^. The disc con¬ 
tains a 5 Mj planet on a circular orbit at a fixed radius of 
rp 40 au. The planet is implemented as a point mass parti¬ 
cle snbject only to the star’s gravity, whereas the SPH gas 
and dust particles feel the gravitational potential of the star 
and planet (see Fouchet et al. 2007). The disc is initially set 
up as in Fouchet et al. (2010) with a constant surface den¬ 
sity Er 19.67 kgm“^ with radial extension from 4 to 120 au, 
a vertically isothermal temperature profile Tr 15r/rp“* K, 
and is free to evolve. Grains have an initially uniform size 
^0 10 yum and are allowed to grow and fragment. In Gonzalez 
et al. (2015), we presented a suite of simulations that vary 
the fragmentation threshold: Vf^g 10, 15, 20, 25 ms“* and oo 
(i.e. growth only), evolved for 100 000 yr or ~400 planetary 
orbits, and discussed our choice of parameters and their in¬ 
fluence on the results. In particular, our disc temperature is 
colder than what observations suggest (e.g. Williams & Best 
2014). Changing its value would affect gap formation (only 
moderately for a high-mass planet), as well as fragmenta¬ 
tion (Vrel X CsX yff). In the midplane, St is independent of 
T, leaving dust dynamics unchanged. In this Letter, we fo¬ 
cus on the Ffjag 15 ms“^ case, for which we see a second 
apparent gap. 

Figure 1 shows the gas and dust volume densities pg and 
Pd, as well as the grain size i in the disc midplane at the end 
of the simulation. The planet at 40 au has opened a gap in 
the gas phase and much of the disc interior to the gap has 
been accreted onto the star, leaving a central low-density 
area. In the dust phase, the gap is more prominent, with 
sharper edges. Some material remains in the inner disc and 
grains have accumulated at the dust trap located at the gap 
outer edge. These properties are commonly seen in studies 
of planet gaps in gas+dust discs (see, e.g., Paardekooper & 
Mellema 2004; Fouchet et al. 2010). An unexpected feature 
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Figure 2. Azimuthally averaged radial profiles of the gas and 
dust volume densities (top), dust-to-gas ratio (middle), and mean 
grain size (bottom) in the midplane after 100 000 yr. The vertical 
dashed line shows the planet orbital radius, the vertical dash- 
dotted lines show the locations of pressure maxima (see Fig. 6) 
and the horizontal dotted line shows the initial dust-to-gas ratio. 



Figure 3. Meridian plane cut of the dust distribution (top) and 
radial grain size distribution (bottom) after 100 000 yr. The color 
represents the volume density and is the same as in Fig. 1. The 
dashed line marks the grain size for which St = 1 in the disc 
midplane (see Sect. 3). 


is the second high-density ring in the dust disc at r~ 90 an 
visible in the density map (centre panel of Fig. 1) and in the 
radial profile of the midplane volume density (Fig. 2). Large 
grains (mm to cm in size) are abundant in the high-density 
regions (right panel of Fig. 1 and bottom panel of Fig. 2), 
where relative velocities between grains are lowest and grain 
growth is more efficient. The dust population is best char¬ 
acterised in Fig. 3, showing the distribution of grains in the 
meridian plane and their radial size distribution. Two popu¬ 
lations can be distinguished: grains that have been trapped 
in a narrow ring at the outer planet gap edge (at r ~ 60 an) 
and have grown, and grains that have started to grow in the 
outer disc, migrated inwards, decoupled from the gas, and 
piled up in a wider ring ~90 an from the star. We refer the 
reader to Gonzalez et al. (2015) for a detailed discussion of 
dust evolution, as well as Sect. 3. 
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Figure 4. Simulated ALMA observations of a disc viewed nearly 
face-on at d 140 pc and S -23° for an integration time of 1 h and 
angular resolution of 0.1". From left to right: A 350 pm, 850 pm, 
and 1.3 mm. The scale on each image is in arcseconds, with the 
beam size represented at its bottom left corner, and the colour bar 
gives the flux in mjy/beam. The planet’s position is marked with 
a cross. 
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Figure 5. Radial cut of the brightness profile in the ALMA syn¬ 
thetic images for A 350 pm, 850 pm, and 1.3 mm. The vertical 
dashed line shows the planet orbital radius and the vertical dash- 
dotted lines show the locations of pressure maxima (see Fig. 6). 

2.2 Synthetic ALMA images 

We produced synthetic ALMA images for the final state 
of our disc. We first computed raw intensity maps from 
the resulting dust distribution using the 3D Monte Carlo 
continuum radiative transfer code MCFOST (Pinte et al. 
2006, 2009). It computes the scattering, absorption and re¬ 
emission of photons by dust grains assumed to be homoge¬ 
neous spheres and composed of astronomical silicates, with 
optical properties derived according to Mie theory. These 
maps were then passed to the CASA simulator for ALMA to 
obtain synthetic images for a given observing conhguration 
(wavelength, angular resolution, integration time). For more 
details about this procedure, see Gonzalez et al. (2012). 

Figure 4 shows the resulting synthetic images at three 
different wavelengths: 350 pm, 850 pm and 1.3 mm. They 
were computed for a nearly face-on orientation, a distance 
d 140 pc and a declination 6 -23° (for which the source 
passes through the zenith at the ALMA site), an integration 
time t 1 h and angular resolution 9 0.1" (found to be optimal 
observing parameters for gap detection by Gonzalez et al. 
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Figure 6. Azimuthally averaged profiles of the gas pressure (top) 
and pressure gradient (bottom) in the midplane after 8, 56 and 
398 planet orbits (~ 2000, 14000 and 100000 yr). The vertical 
dashed line shows the planet orbital radius and the vertical dash- 
dotted lines show the locations of pressure maxima at the end of 
the simulation. The inset shows a larger vertical scale. 

2012). We also plot in Fig. 5 radial cuts of the brightness 
profile for the three wavelengths. All images show the central 
emission from the dust in the inner disc (interior to 20 an, 
see Fig. 1), the deep and wide planet gap (centred on 40 au), 
a bright ring corresponding to the outer edge of the planet 
gap (near 60 au), a second shallower and narrower apparent 
gap (centred on ~ 75 au, which happens to be near the 2:5 
resonance with the planet), and finally a second fainter ring 
(centred on ~ 90 au). The synthetic observations thus recover 
all features of the dust distribution. Note that a higher disc 
temperature would increase the computed fluxes, but would 
preserve the gap structures (see Sect. 2.1 for the influence 
of T on the dynamics). 


3 DISCUSSION 

We hnd that, in a CTTS disc containing a 5 Mj planet 
on a fixed circular 40 au orbit, a fragmentation threshold 
of Ffrag 15 ms“* produces two almost axisymmetric gap 
and ring pairs in the corresponding synthetic images at 
(sub)millimetre wavelengths (Fig. 4). What is their origin? 
Lower continuum emission at a given location corresponds 
to lower dust opacity at millimetre wavelengths, which can 
be caused by a local minimum of the dust density or by grain 
sizes outside of the range contributing significantly at these 
wavelengths (i.e. 100 ^m to a few mm, for the dust optical 
properties assumed here), or both. It is important to note 
that the two gaps in our disc have different origins. 

The planet is massive enough to carve a gap in the gas 
phase and produce a local pressure maximum at the outer 
gap edge (r ~ 60 au) very early in the disc evolution - see the 
midplane pressure and pressure gradient prohles in Fig. 6. 
Dust particles are trapped at the pressure maximum, accu¬ 
mulate, and grow more efficiently as they concentrate. This 
so-called particle trap generates the first ring just outside 
the planet gap (see also Pinilla et al. 2012). 

The second ring, on the other hand, is caused by a grain 


pile up rather than a trap due to the planet. Grains can grow 
when Vj-ei (which scales as oc in our disc) is small 

enough to stay below Ffrag, i-S- in the disc outer regions. 
The growth timescale varies as (Laibe et al. 2008), it 
is thus shorter at the inner boundary of this growth region. 
Growing grains migrate inwards, more rapidly as their size 
approaches the value for which St 1 (plotted as a dashed line 
in Fig. 3 at the end of the simulation, see Gonzalez et al. 
2015 for the time evolution). As grains continue to grow, 
they overcome the St 1 regime, causing them to decouple 
from the gas, which slows their radial migration and results 
in a pile up near 90 au. This enhances the local concentra¬ 
tion of grains, and the dust-to-gas ratio approaches unity 
(see the middle panel of Fig. 2), which i) stops grain mi¬ 
gration in the accumulation region (Nakagawa et al. 1986) 
and ii) starts to modify the gas density profile due to the 
backreaction of dust on gas. Indeed, the drag force between 
gas and dust is symmetric and in a dust concentration, it 
is gas that is dragged by the dust. This effect can only be 
seen in self-consistent simulations of the simultaneous evo¬ 
lution of gas and dust taking backreaction into account. In 
other words, the complex interplay between growth and mi¬ 
gration induces dust self-accumulation, which causes the gas 
to accumulate at the dust concentration and form a second 
gas pressure maximum around r ~ 90 au at late times (see 
Fig. 6). This effectively results in a second particle trap, 
causing the second ring. This mechanism, predicted theo¬ 
retically by Laibe (2014), is independent of the presence of 
a planet, and was also seen in the very inner disc for the 
growth-only simulations in Gonzalez et al. (2015). Here, we 
identify for the first time a self-induced dust pile up in sim¬ 
ulations taking fragmentation into account. 

Counter-intuitively, the value of the fragmentation 
threshold plays an essential role in the generation of this 
second ring. The dust pile up occurs near the location 
where St 1 for the largest grains, for which Vrel is also 
equal to Vfrag, i.e. where Csy/a^ ~ For a sound speed 

given by Cs Cs,0^/^0 this location is approximated by 
rpu/ro ~ In our disc and for Vfrag 15 ms^', 

^pii ~ 90 au. For larger values of Vfrag, grains would pile up at 
smaller orbital distances (~ 50 and ~ 30 au for Vfrag 20 and 
25 ms^*) but they are trapped at the planet gap outer edge 
before they can decouple. There is only one ring in this case. 
On the other hand, for smaller Ffrag, only grains at very large 
orbital distances can grow, but because of their even lower 
growth rate they do not reach St 1 and do not decouple be¬ 
fore the end of the simulation (see Gonzalez et al. 2015, for 
illustrations of both cases). Vfrag ranges from 1 and 10 ms”* 
for compact silicate and icy grains (Blum & Wurm 2008) 
to several tens of ms”* for porous aggregates (Yamamoto 
et al. 2014). The value of 15 ms”*, producing a second ring 
that is well separated from the outer gap edge, would match 
moderately porous ice. However, different values of Vfrag, or 
similarly of Uss or the disc temperature, can still lead to two 
separate rings if the planet location is changed. 

The second apparent gap is therefore the region of lower 
dust density between both dust accumulations. Its proper¬ 
ties differ from those of a planet gap: the decrease in the 
density of both gas and dust is modest and the grains have 
grown to moderate sizes. However, the observational signa¬ 
tures of both the planet gap and the apparent gap are very 


MNRAS 000, 1-6 (2015) 




































ALMA images of discs: are all gaps carved by planets? L5 


similar, both in the images (Fig. 4) and in the brightness 
profiles (Fig. 5), so that the shallower apparent gap seems 
to mimic a gap that would be carved by a lighter planet 
orbiting at ~ 75 an (all the more so that it would lie close 
to the 2:5 resonance with the planet at 40 au). Measuring 
this gap’s width would give an estimate of this hypothetical 
planet’s mass (Crida et al. 2006), but would not unambigu¬ 
ously prove its existence. Cuello (in preparation) ran a series 
of simulations with two planets embedded in the same disc 
as ours. With planets of 5 Mj at 40 au and 1 Mj at 80 au, the 
outer planet gap has a similar depth In the gas as our ap¬ 
parent gap (also in agreement with Kanagawa et al. 2015) 
but a much deeper gap in the dust. Conversely, a lighter 
outer planet would be needed to reproduce the depth of our 
apparent gap in the dust, but it would carve an even shal¬ 
lower gap in the gas. If a suitable tracer of gas In the gap 
can be found, ALMA observations of both gas and dust and 
the comparison of gap depth in both phases would thus be 
able to distinguish an apparent gap from a planet gap. 


4 CONCLUSION 

We have presented 3D hydrodynamical simulations of the 
evolution of gas and dust in a CTTS disc containing a 
planet which carves a gap in the disc. We self-consistently 
treat the interaction of gas and dust via aerodynamic drag 
and include grain growth and fragmentation. We then com¬ 
puted synthetic ALMA images from the resulting disc struc¬ 
ture. We have shown that when the fragmentation threshold 
Vfmg 15 ms'b the dust dynamics results in the first docu¬ 
mented self-induced pile up of grains when fragmentation is 
included, at large distances from the star. This forms a dense 
dust ring containing mm-sized grains in the outer disc, away 
from the planet gap and the dust trap at its outer edge. In 
the images, in addition to the gap+ring pair corresponding 
to the planet gap and its dense outer edge, this feature pro¬ 
duces a second ring separated from the outer gap edge by 
an apparent gap, effectively resulting in a second gap-fring 
pair. This second, exterior gap does not contain a planet. 

The recent spectacular ALMA image of the HL Tau disc 
(ALMA Partnership et al. 2015) shows as many as seven gap 
and ring pairs. Is each of them an actual planet gap? The 
answer is probably no and we have shown here one example 
for which it is not the case. More gaps in discs are likely to 
be detected as ALMA is ramping up to full operation and 
they will have to be interpreted with care. Comparison of 
gap depths in both gas and dust would help to distinguish 
a planet gap from an apparent gap such as presented here. 
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